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Abstract 

The problem of disturbed flow in the laminar boundary layer is in- 
vestigated when a heating element is located on the surface of the body. 
The flow is supposed to be laminar one. It is shown that the problem 
may be solved in terms of free interaction theory. The solution of linear 
flat problem is constructed. The results of asymptotic analysis are also 
presented. 



1 Introduction 

The problems of boundary layer control remain important being connected 
with e.g., the increase of fuel prices. The reduction of vehicle resistance ac- 
complished in particular by means of this control arrives at the reduction of 
fuel consumption and thus, to the increase of the flight effectiveness. Among 
a number of methods for boundary layer control the method in which heating 
elements are located on the surface and vary the temperature in the boundary 
layer, plays the crucial role. The use of these elements is similar, in some 
sense, to the application of small humps on the surface which also influences 
the flow in the boundary layer. In this case the problem of influence of local 
heating of the surface on the flow in the boundary layer appears. The crucial 
characteristics in this problem are pressure and shear-stress behaviour. The 
perturbation of these function in the neighbourhood of the heating element 
determines how to influence the boundary layer to delay laminar-turbulent 
transition, reduce the resistance etc. The general consideration of the problem 
implies the construction of solutions of the nonlinear boundary problems and 
investigation of perturbations depending on time. However, the approach on 
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the current stage requires to consider simpler stationary models and can be 
linearized. This approach enables to state the general mechanisms of the flow 
adjacent to the heated region and estimate the physical effects of action of 
such heating on the flow. 



2 Scales evaluation and problem statement for 
the viscous sublayer 

We consider the subsonic or supersonic^ compressible flow over the semi 
infinite flat plate when the flow still may be considered laminar but Remolds 
number is sufficiently large. Denote Re = e~ 2 . Thus, e is assumed to be small 
parameter. In general, the flow is supposed to be two dimensional and nonsta- 
tionary one. Preliminary, let us transform the original Navier-Stokes equations 
to nondimensional ones. Introduce the appropriate nondimensinal variables 
denoted with the same letters xl,yl, zl,lu^t, u^u, u^v , PooP, v^R"^, PooP) 
poouioP- Here Uoo, Poo, Poo are a longitudinal velocity, density and the coefficient 
of viscosity in the undisturbed flow over the surface where Euler equations are 
valid, I, distance from the leading edge to the place on the plate where the 
heating region is located. Substituting these functions in Navier-Stokes equa- 
tions, we obtain the nondimensional equations differing from the former ones 
only by the multiplier I /Re in the right sides of the equations. 

In the problem under consideration there are three different regions with the 
specific equations. These regions naturally appear while asymptotic expanding 
of solutions of the Navier-Stokes equations and balancing the orders of the 
convective, dissipative and pressure terms. This structure is usually called 
triple-deck. The viscous sublayer on the bottom of the boundary layer is called 
inner deck, the boundary layer itself and the inviscid flow are called the main 
deck and the outer deck respectively. We will use these terms interchangeably 
with the terms the viscous sublayer, boundary layer and inviscid flow. 

In previous works the problem of temperature and catalicity jumps has been 
systematically investigated [U [2] • In this case we assume that the temperature 
disturbances appear in a region with finite length to be estimated. It was shown 
in[U [2] that in the problem with the temperature jump AT ~ T ~ 0(1) and 
the disturbances propagate on the distance Ax ~ e 3//4 . Therefore we assume 
that the order of the length of the heating region is the same as that of the 
region of disturbances propagation i.e., a ~ e 3//4 where a is the length of heated 

1 We consider the subsonic problem but the discussion in the current section is general 
both for subsonic and for supersonic case. 
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region. It was shown in [3] that this situation leads to forming of a local 
effective irregularity (hump) on the surface and the problem becomes similar 
to that of compressible flow over the local irregularitiesjUElElITj. Consider the 
region where the convective, dissipative, and pressure terms of Navier-Stokes 
equations are of the same order. The estimations of the transversal size of the 
interaction region follow from the asymptotic equivalence of the convective and 
dissipative terms 

du 2 d 2u 
dx dy 2 

Taking into account that in the sublayer (inner deck) u ~ y/e we obtain the 
estimate y ~ e 5 / 4 . Similarly, the asymptotic relation 

dp 2 d 2 u 
dx dy 2 

leads to the estimate of the pressure in the inner deck p ~ e 1 / 2 . Time changing 
scale is represented by r ~ e^ 1 ^ 2 . 

The above estimates show that under these circumstances free interaction 
mode is valid [HI El CD]- Asymptotic representations in the regions located above 
the inner deck were derived in [3]. Consider in detail the boundary problem in 
the viscous asymptotic sublayer of the nonlinear disturbances which is located 
beneath the asymptotic Prandtl boundary layer and outer inviscid flow. The 
functions in the sublayer denoted by the index 3 have the following asymptotic 
expansions 

x = 1 + ax 3 , y = a l ^ty 3 , t = rt 3 

u = a l ^u 3 + . . . , v = ea~ 1 / 3 v 3 + . . . , 

p = Poo + a 2/3 p 3 + ..., T = T 3 + ..., p = p 3 + ..., 

where a is the longitudinal scale estimate for the heating region. Substitut- 
ing the above expansions in Navier-Stokes equations and letting e —>■ we 
arrive at the set of equations for the viscous nonstationaty layer which is two 
dimensional nonstationary boundary layer 

Qi dp 3 dp 3 u 3 dp 3 v 3 

(sh— + u — +v — ^ + — = — ( — ^ 

P3 \ dt 3 Us dx 3 V3 dy 3 J dx 3 dy 3 \dy 3 ) 

(2.1) 

3 V dt 3 3 dx 3 3 dy 3 ) dy 3 \Pidy 3 ) 
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— = 0, p 3 T 3 = 1, 

In these equations Pr is Prandtl number, Sh = r~ 1 e 1 / 2 , Struhal number. The 
boundary conditions are formulated on the basis of no slip conditions and 
matching the solutions with these in other regionsplj. For the pressure the 
interaction condition holds having the following form for the subsonic case [3] 

+00 

B 1 r dD/di 



P: ' = - WT^W J T^*' 

—00 

where D is an arbitrary function. The initial conditions are derived by match- 
ing with the wall adjacent part of the undisturbed boundary layer on the plate 

u 3 -> Ay 3 ,T 3 -> T 20 (y 3 ),p 3 ,D -> 

as x 3 — > — 00. On the surface of the plate y 3 — and the no slip conditions 
yield u 3 = v 3 = 0,T 3 = T 3w . At a large distance from the plate y 3 — > 00 and 
the boundary conditions have the form u 3 — > A(y 3 + D),T 3 — > T 20 (0), where 
T 2 o is the main term of temperature asymptotic expansion in the main deck. 
Let us introduce the additional transform of the independent and dependent 
variables to simplify the numerical analysis of the problem 

4-1/8 r *»> 

x 3 = x b , y 3 = A 1 / — , 

Jo Pb 

u 3 = A 2 ^ 3 u b , 

P3 = P2o(0)A 2 / 3 p b , P 3 = p 20 (0)p fe , T 3 = T 20 (0)T b 

B 2 = fi 1 7r- 1 p2o 1 (0)^ _5/3 (l - M 2 )- 1/2 , D = A~ 1/3 d 

The value v b is established depending on the new variables. In addition, set 
Pr = const and p 3/ u = 1 implying that the coefficient of viscosity is linearly 
dependent on the temperature. Thus, the boundary problem (12. ID is rewritten 
as follows 

dub | dvb _ q 
dx b dy b 



du b du b du b dp b d 2 u b 
c% dx b dy b dx b dyf 



i> 

(2.2) 
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dT b dT b dT b d 2 T b 
dt b dx b dy b dy£ 
with the boundary conditions 



u b (x b , 0, t b ) = v b (x b , 0, t b ) = 0, 
T b (x b ,0,t b ) = T w (x b ,t b ), 

P+OO 

u b -> Ub + d, di = / (1 - T b )drj + d, T b (x, oo) — >• 1, y b -> oo 
Jo 



+oo 



d(-oo) = 0, Pb(x b ,t b ) = -B 2 ! ^p—^- dC, 

J Q-Xb 

— oo 

Note that the interaction condition can be transformed into the other form 
as[5j 

J X b q 



where 

- + 0O 



r+oo 

Aib(x b ) = (1 - T b )dr] + A b , u b ^y b + A b (x b ),y b -> oo. 
The last relation we will use in the following sections. 



3 Construction of the solution of flat linearized 
stationary problem 

3.1 General construction of the solution for the case of 
arbitrary temperature in the viscous sublayer 



Let us consider the stationary problem. Then (12.2ft may be transformed as 
follows 

du b dvh _ g 
dx b dy b 



du b du b dp b d 2 u b 
u b T. h v b - h T b — = — -=- (3.1) 

dx b dy b dx b dyl 



dT b dT b d 2 T h 



dx b dy b dy% 

Assume that the temperature disturbances are small i.e., AT << T ~ 0(1). 
Then the solution of (13.11) may be constructed by expanding the flow functions 
in the main undisturbed part of the boundary layer. Let us introduce the small 
parameter A = AT and represent the flow functions as follows 

u b = y b + XU + 0(X 2 ), v b = XV + 0(X 2 ), A lb = XA l + 0{X 2 ) 

T b = 1 + AT + 0(A 2 ), p b = XP + 0(X 2 ), A b = XA + 0(X 2 ) 

as A — > 0. Substitution of these expansions in (13. ip and boundary condition 
and equalization of the terms with the same powers of A yield the set 

dU dV 
dx b dy b 

dU dP d 2 U 

dx b dx b dyl 

dT d 2 T 



Vb- 



'dx b dyl 
The boundary conditions take the form 

U{x b ,0) = V(x b ,Q) = 0,T{x b ,0) = T w (x b ), 

r+oo 

U^A(x b ), Ai(x b ) = Tdrj + A(x b ), T(x b , oo) -> 0, y b -> oo 

J 

+oo 

P'(0 



A(-oo) = 0, A'l{x b ) = - J 



x b -i 

— oo 

Applying Fourier transform with respect to x b to the last set, we obtain 

i£U + v' = o, 

i£y b U + V + i£P = U", (3.2) 
i£y b f = f" 

Consider the last equation of (13. 2p . Set Y = (i^Y^ 3 y b and y b > 0. Then the 
equation becomes 

f" = YT, 
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where primes denote differentiating with respect to Y. Thus, in the result 
of variable changing the variable Y is complex and the equation for Fourier 
image is Airy equation. To determine the values of Y uniquely it is necessary to 
specify the branch of multi- valued function (z£ ) 1//3 . Let us consider the complex 
plane to = 77 + i£. We find for Airy function Aif^u; 1 / 3 ^). Cut the plane to along 
the negative part of real half axis. In the plane with the cut the function to 1 / 3 
allows to specify regular branches. Then take the branch of cubic root such 
that at uj = 1 we have u 1 ^ 3 = 1. In the following consideration we will use the 
consequence of previous procedure, that for the branch under consideration 
Ai((i0 1/3 2/&) = Ai(^/ 3 y b e iw / 6 ) as f > and Ai((^) 1/3 2/&) = Ai^ 1 ^ e~ i7r / 6 ) as 

£<o. 

The boundary condition at infinity is derived from those for (I2.2p and has 
the form T(oo) = 0. The second boundary condition is set using the temper- 
ature distribution on the wall and taking into account its Fourier transform 
and has the form T(£) = /(£). Using the new variable Y the equation for 
Fourier images of temperature becomes Airy equation. The boundary condi- 
tion at infinity shows the solution of Airy equation must be chosen recessive at 
infinity. This solution depends on Y complex plane segment under considera- 
tion. However, the above suggested specification of the branch shows that in 
the case under consideration the argument of Airy functions is located in the 
segment defined by |(| argY) < n/3 where the recessive solution is represented 
by means of Ai(Y) and hence the solution is given by 

f(Z,y b ) = K(£)M(Y). 

The boundary condition in zero point gives 

ho 



Ai(0) 



and we obtain the solution for the images of the temperature. 

Consider two other equations of (13.21) . Differentiating the second one with 
respect to y& and substituting the result in the first equation we obtain 

U'" = -i£y b U'. 

Introducing notation U' = F and changing variables as above, i.e. Y = 
(iCY^Vb, we arrive at Airy equation 

F" = YF. 

In the last equation the primes imply differentiation with respect to Y. The 
boundary condition for this equation follows from the asymptotic relation for 
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the longitudinal velocity as y b — > +00. The solution of the equation satisfying 
the boundary condition is represented by Airy function Ai(Y). The solution 
also depends on the complex plane segments as in the case of the temperature 
equation and we do not repeat the above analysis. Notice that the velocity 
distributions are not crucial functions in the problem under consideration. It 
is necessary to obtain the wall-shear and the pressure in the viscous sublayer. 
Thus, we find 

= M(£)M(Y). (3.3) 

To determine M(£) we observe that the longitudinal impulse equation 
yields at y b — 

dx b dyl %b) 

For the pressure derivative we denote P' x = G(xf ) ). Then the equation for 
images takes the form 

G = U"(£,0) (3.4) 

Differentiating (13.31) with respect to y b , evaluating it at y b — and substituting 
in (13.41) . we have 

G = M(O(*O V3 Ai'(0). 

The second necessary relation for evaluating G and M(£) follows from the 
interaction condition and has the form 



+00 +00 

P'(C) 



U"(x b , 00) + J T"dr] = - J 



x b -C 



d( 



Applying the Fourier transform to the last equality, we obtain 

+00 +00 +00 

(z£) 2 ^(£, oo) + (*0 2 J T(e, V)dv = - J e"^ dx b J 



+00 



. i7r M J e - i ttp'(()d( = iir^G 



—00 

Next, observe that from the foregoing analysis follows 



(i0 1/3 J Ai(0) 
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we arrive at the following condition 





where we use the notation k = (- Ai'(O)/ J +o ° ki(z)dzf /A = 0.827[12j. Here it 
is necessary to return to the problem of Airy function argument choosing. The 
last discussion arrive at improper integral of Ai(y) calculation which converges 
as |arg(y)| < 7r/3[TT]. The branch of u 1 ^ 3 must be specified to provide the 
convergence of the integral both as £ > and £ < 0. The above specified 
branch provides this convergence. 

Taking into account two conditions for G and M(£) we obtain 

d= Ai'(0) (i0 5/3 /(0 



Ai(0) ^1^4/3 + (^)4/3 
Ai(0) C i7r Mjfe4/3 + )4/3 



M(0 

,-U(V) ( i7T 

Lastly, for the Fourier images of the shear we find from the previous relation 
r = U' yb (d, 0) = M(0 Ai(0) = — |f f° / 4/3/(0 ; 

In all above formulas the plus is chosen at £ > and the minus at £ < 0. 

3.2 Construction of the solution and asymptotic 
estimations of flow functions 



The relations for the functions M(£) and G show that in the case under 
consideration the influence of the function /(£) on convergence of integrals 
for the wall-shear and the pressure is more significant than in the supersonic 
case[l3]. It follows that the shape of local heating region has more influence 
on the functions in the lower-deck than in the supersonic case. Let us note 
the properties of /(£) to be used below. Since we consider the local regions 
of heating, we assume this function to be Fourier transform of some finite 
function. Then, the original finite function is assumed to be even, implying 
its Fourier transform to be a real function. We consider the pressure case in 
more detail for the below analysis is similar for the wall-shear. 
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The pressure gradient on the wall is given by 



+00 



P' 



G(x b ) 



2tt 



1 Ai'(O) 
"2n Ai(0) 



(<e) 5/s /(0 



+00 



-ivrfc 4 / 3 + («0 4/3 



<r : 7(0 



Accomplishing the simple transformations and using the results of §3.1, intro- 
ducing the variable r > such that as £ > we have i£ = r e l7T / 2 and as £ < 
= re"™, we obtain 



G(x b ) 



^Ai^O) 
'2tt Ai(0) 



-i5n/6 , 



,5/3 /( . 



-ink 4 / 3 + e ~ i27T / 3 r 4 / 3 



-irxi, 



dr + 



+ 



a i57r/6 , 



■5/3 J(, 



ink 4 / 3 + e i2,r / 3 r 4 / 3 



(3.5) 



The heating region of the surface has some length, therefore to estimate the 
shear and the pressure we set the temperature on the wall by means of the 
following finite even function 



T(x b ,0) 



0.2, |ar 6 | < 0.5 
0, \x b \ > 0.5. 



Fourier transform of this function yields 

/(0 = T sm (- 



It follows that /(£) — > as £ — > 00 and even. In addition, it is easily seen 
that in (13.51) we obtain the sum of two conjugate functions. Substituting the 
last relation in ( 13. 5ft and taking into account the complex conjugation, we 
obtain the integral representation of the general solution in the sublayer for 
the pressure 



P'(x b ) = - 



0.4Ai'(0) f r 2 / 3 sin(rx b — ip) sin(r/2)dr 

7T Ai(0) J v /^8/3 + y3 7rfc 4/3 r 4/3 + r 8/3 ' 



(3.6) 
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Figure 1: Pressure in the viscous sublayer. Dashed lines - asymptotic curves 



where ip is determined by 

/ irk 4 ' 3 + v/3r 4 / 3 ' 

w = arctan — ■ = - 

\r 4 /3 + ^A; 4 / 3 

In the similar way for the shear we obtain 

+00 

dU 0.4 f r 1 ' 3 sinfrxf, — 0) sin(r /2)dr 
— (x b ,0) = 



8y b 7T J v / 7r 2/ |; 8/3 + v /3 7rfc 4/3 r 4/3 + r 8/3 ' 





where 

6 = arctan 



The results of pressure estimation are depicted in fig. 1, and the wall-shear, 
in fig. 2. The pressure estimation is obtained by means of integration of (13. 6h . 

Let us consider the asymptotic behaviour of flow functions and firstly con- 
sider the wall-shear case. The main relation we represent as the sum of integrals 
as follows 

qjj q 4 + Tr 1 / 3 (^-iTk^ 3 + r 4 / 3 j sin(r/2) cos(rx b )dr 



d V ~S Xh ^ TT J ^2^8/3 + ^3^4/3^/3 + r 8/3 + 
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Figure 2: Wall-shear in the viscous sublayer. Dashed lines - asymptotic curves 



+oo 



0.4 f r 1 / 3 !^ 4 / 3 sin(r/2) sm(rxb)dr 



0.4 1 

+oo 



vr J tt 2 A; 8 / 3 + ^37rfc 4 / 3 r 4 /3 + r 8/3 



+oo 



+ OO +oo 

J r 1 / 3 h 1 (r)e ir ( Xb+1 Wdr+ J r^h^r) e -^(.x b -i/2) dr _ 



+■ 



0.4 1 

7T 4 



+oo +oo 

J r^hir) e^-W dr - J r^h^r) e ~ ir ^ +1 ^ dr 



+oo +oo 

J r^ 3 h 2 {r) e -(^ +1 /2) dr _ J r ^ h ^ r ) e -<r(x 6 -i/2) dr _ 



+oo 



+oo 



- J r 1/3 h 2 (r)e ir{xb - 1/2) dr+ J r 1/3 h 2 (r) e - ir ^ +1/2) dr 



where we use the following notation 

yJttA.4/3 _l r 4/3 



Ink 4 / 3 



7T 2 /fc 8 / 3 + VSnk 4 / 3 r 4 / 3 + r 8 / 3 
Observe that 



7i 2 k s / 3 + \/37rA; 4 / 3 r 4 / 3 + r 8 / 3 ' 



^(0) = ^vrA; 4/3 , h 2 (0) = ]-7ik 4/3 , h[(0) = h' 2 (0) = 0. 
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For the above integrals the asymptotic relations as Xb — -> +00 can be given 
being the specific case of relations described in [TJ], namely 

/ r 1/3 h(r) e ±ir5 dr ~ ^ #- fc -i(f , 0)/i (fe) (0), £ -> +00 
fc=0 

where 

o) = ^^H)r (* + r fc "*, 

and represent the integrated kernels of the main kernel K (^,r) = r 1 ^ 3 e"'^. 
Even though the main kernels Xo(6, r ) are different for each of the integrals but 
being of oscillatory type they allow to construct the asymptotic expansion [14j. 
Applying the suggested expansion to each of the above integrals, we obtain 
the following estimation of the downstream wall-shear 

As Xb — > —00 the slight modification of the analysis yields the following asymp- 
totic relation 

Tfcl** °> = -3^F7s r (5) (Z^m + o«-^)- ,/s )). * - — ■ 

Thus, the wall-shear asymptotics upstream and downstream has the same 
order of vanishing differing this case from the supersonic one. The obtained 
asymptotic curves are depicted on the fig. 2, and fit well with the precise 
solution at large Xb- Note that in the above general asymptotic formula the 
function h(r) is supposed to be infinite differentiable in zero, which can cause 
difficulties as e.g., in the case under consideration. However, the main formula 
relies on the integration by parts and therefore may be applied successively till 
the smoothness of h(r) allows. Nevertheless, we also note that the smoothness 
conditions on h(r) in zero may be significantly weakened, not affecting the 
obtained results. 

The asymptotics of the pressure is estimated quite similarly to the case 
involved and we represent only the final results. Estimating in the similar way 
for ( 13. 6p we obtain 

[Xb) 9vr 2 ^/3 Ai(0) x 8/ 3 + ) ' b + 
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Since the integration of this asymptotic expansion is valid[TT], we obtain the 
estimate 

„, N 2r(|) Ai'(O) 1 / 1 \ 

P{Xb) = iSP Ai(0) ^ + ° U^J ' Xb ^ + °°- 

The similar consideration for the case Xf, —>■ — oo gives the estimate 

p, Xb) F (l) A '(°> 1 t: ( 1 ) z^-oc 

On the fig. 1 the curves of asymptotics and the precise solutions of the pressure 
are depicted. 



4 Discussion 

The results show that though qualitatively the flow functions are similar 
to those for the problems of local irregularities but there are also significant 
differences. The curves shown in fig. 1,2 allow to conclude that the wall-shear 
and pressure are described by continuous functions. In addition, the shear in 
the problem under consideration is similar to that obtained in the problem 
for the supersonic flow|13j. Downstream the shear increases(fig. 2) implying 
favourable pressure gradient (fig. 1). The possible points of flow separation are 
boundaries of heating region upstream and downstream. However, the asymp- 
totic analysis reveals the vanishing of perturbations in the problem is different 
from that in the supersonic flow. Particularly, for the wall-shear in the su- 
personic case it was obtained that the perturbations upstream vanish stronger 
than in the subsonic case. This fact is easily explained from the physical point 
of view. The supersonic flow "blows off" the upstream perturbations while in 
the supersonic case the flow allows the perturbations to propagate upstream. 
Mathematically this fact can be explained by pole of the integrand for the 
flow functions in the supersonic case in the lower half plane resulting in ex- 
ponential vanishing upstream. The similar situation is valid for the pressure. 
As we mentioned above, the shape of the effective hump in the problem un- 
der consideration strongly influences on the pressure in the subsonic case that 
also can be explained by the fact that the perturbations propagate further up- 
stream. Asymptotic analysis also shows that the obtained representations for 
the asymptotic behaviour of flow functions fit well with the precise formulas on 
the distance to be equal to about the half of this region and can be used in the 
numerical estimations e.g. for the numerical approximation of the interaction 
condition. 
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